非线性玻尔兹曼方程的傅里叶谱方法

玻尔兹曼方程作为空气动理学中最基本的方程之一,是连接微观牛顿力学和宏观连续介质力学的重要桥梁.该方程描述了一个由大量粒子组成的复杂系统的非平衡态时间演化:除了基本的输运项,其最重要的特性是粒子间的相互碰撞由一个高维,非局部且非线性的积分算子来描述,从而给玻尔兹曼方程的数值求解带来非常大的挑战.在过去的二十年间,基于傅里叶级数的谱方法成为了数值求解玻尔兹曼方程的一种很受欢迎且有效的确定性算法.这主要归功于谱方法的高精度及它可以被快速傅里叶变换加速的特质.本文将回顾玻尔兹曼方程的傅里叶谱方法,具体包括方法的导出,稳定性和收敛性分析,快速算法,以及在一大类基于碰撞的空气动理学方程中的推广.     

NIFTy 3 – Numerical Information Field Theory: A Python Framework for Multicomponent Signal Inference on HPC Clusters

NIFTy , “Numerical Information Field Theory,” is a software framework designed to ease the development and implementation of field inference algorithms. Field equations are formulated independently of the underlying spatial geometry allowing the user to focus on the algorithmic design. Under the hood, NIFTy ensures that the discretization of the implemented equations is consistent. This enables the user to prototype an algorithm rapidly in 1D and then apply it to high‐dimensional real‐world problems. This paper introduces NIFTy  3, a major upgrade to the original NIFTy  framework. NIFTy  3 allows the user to run inference algorithms on massively parallel high performance computing clusters without changing the implementation of the field equations. It supports n‐dimensional Cartesian spaces, spherical spaces, power spaces, and product spaces as well as transforms to their harmonic counterparts. Furthermore, NIFTy  3 is able to handle non‐scalar fields, such as vector or tensor fields. The functionality and performance of the software package is demonstrated with example code, which implements a mock inference inspired by a real‐world algorithm from the realm of information field theory. NIFTy  3 is open‐source software available under the GNU General Public License v3 (GPL‐3) at https://gitlab.mpcdf.mpg.de/ift/NIFTy/tree/NIFTy_3 .     

On the Relative Gain Array (RGA) with singular and rectangular matrices

This paper identifies a significant deficiency in the literature on the application of the Relative Gain Array (RGA) formalism in the case of singular matrices. Specifically, it is shown that the conventional use of the Moore–Penrose pseudoinverse is inappropriate because it fails to preserve critical properties that can be assumed in the nonsingular case. It is then shown that such properties can be rigorously preserved using an alternative generalized matrix inverse.     

Enhancing signal-to-noise ratio and resolution in low-field NMR relaxation measurements using post-acquisition digital filters

The traditional way to enhance signal-to-noise ratio (SNR) of nuclear magnetic resonance (NMR) signals is to increase the number of scans. However, this procedure increases the measuring time that can be prohibitive for some applications. Therefore, we have tested the use of several post-acquisition digital filters to enhance SNR up to one order of magnitude in time domain NMR (TD-NMR) relaxation measurements. The procedures were studied using continuous wave free precession (CWFP-T₁) signals, acquired with very low flip angles that contain six times more noise than the Carr–Purcell–Meiboom–Gill (CPMG) signal of the same sample and experimental time. Linear (LI) and logarithmic (LO) data compression, low-pass infinity impulse response (LP), Savitzky–Golay (SG), and wavelet transform (WA) post-acquisition filters enhanced the SNR of the CWFP-T₁ signals by at least six times. The best filters were LO, SG, and WA that have high enhancement in SNR without significant distortions in the ILT relaxation distribution data. Therefore, it was demonstrated that these post-acquisition digital filters could be a useful way to denoise CWFP-T₁, as well as CPMG noisy signals, and consequently reducing the experimental time. It was also demonstrated that filtered CWFP-T₁ method has the potential to be a rapid and nondestructive method to measure fat content in beef and certainly in other meat samples.     

Two weight commutators in the Dirichlet and Neumann Laplacian settings

In this paper we establish the characterization of the weighted BMO via two weight commutators in the settings of the Neumann Laplacian ${\mathrm{\Delta }}_{N_{+}}$ on the upper half space ${\mathbb{R}}_{+}^n$ and the reflection Neumann Laplacian ${\mathrm{\Delta }}_N$ on ${\mathbb{R}}^n$ with respect to the weights associated to ${\mathrm{\Delta }}_{N_{+}}$ and ${\mathrm{\Delta }}_N$ respectively. This in turn yields a weak factorization for the corresponding weighted Hardy spaces, where in particular, the weighted class associated to ${\mathrm{\Delta }}_N$ is strictly larger than the Muckenhoupt weighted class and contains non-doubling weights. In our study, we also make contributions to the classical Muckenhoupt–Wheeden weighted Hardy space (BMO space respectively) by showing that it can be characterized via the area function (Carleson measure respectively) involving the semigroup generated by the Laplacian on ${\mathbb{R}}^n$ and that the duality of these weighted Hardy and BMO spaces holds for Muckenhoupt $A^p$ weights with $p\in (1,2]$ while the previously known related results cover only $p\in (1,\frac{n+1}{n}]$. We also point out that this two weight commutator theorem might not be true in the setting of general operators L, and in particular we show that it is not true when L is the Dirichlet Laplacian ${\mathrm{\Delta }}_{D_{+}}$ on ${\mathbb{R}}_{+}^n$.     

Lamb波频域逆散射全聚焦成像方法研究

Lamb波因其传播距离远、衰减小常被用于板状结构的无损检测中,在基于Lamb波损伤检测的诸多成像技术中,全聚焦方法(Total Focus Method,TFM)方法因其成像分辨率高、信噪比高而受青睐。然而Lamb波的频散效应导致时域延时量不能被准确计算,进而影响传统TFM方法对损伤定位及成像的精度;此外,既有的TFM方法仅以回波幅值作为成像指标,忽略了Lamb波与损伤的相互作用,故而不能通过损伤表面的物理参数增强成像质量。针对这两个问题,本文首先在时域TFM基础上发展了频域TFM,在计算中纳入频散关系以规避频散的影响;其次以包含明确的损伤特征参数--反射率为成像指标,结合频域TFM方法建立损伤逆散射模型,以实现对损伤的准确成像。仿真和实验结果表明:频域逆散射TFM成像方法能够实现对铝板结构中的损伤检测,在工程实践中具有应用价值。     

Nonlinear dynamical magnetosonic wave interactions and collisions in magnetized plasma

The one-dimensional nonlinear dynamical wave interactions in a system of quasineutral two-fluid plasma in a constant magnetic field are investigated.The existence of the travelling wave solutions is discussed.The modulation stability of linear waves and the modulation instability of weakly nonlinear waves are presented.Both suggest that the Korteweg-de Vries(KdV) system is modulationally stable.Besides,the wave interactions including the periodic wave interaction and the solitary wave interaction are captured and presented.It is shown that these interacting waves alternately exchange their energy during propagation.The Fourier spectrum analysis is used to depict the energy transformation between the primary and harmonic waves.It is known that the wave interactions in magnetized plasma play an important role in various processes of heating and energy transportation in space and astrophysical plasma.However,few researchers have considered such magnetohydrodynamic(MHD) wave interactions in plasma.It is expected that this work can provide additional insight into understanding of behaviors of MHD wave interactions.